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ABSTRACT 

The X-ray spectrum of the Seyfert 1.5 source NGC 6860 is among the most 
complex of the sources detected in the Swift Burst Alert Telescope all-sky sur- 
vey. A short XMM-Newton follow-up observation of the source revealed a flat 
spectrum both above and below 2 keV. To uncover the complexity of the source, 
in this paper we analyze both a 40 ks Suzaku and a 100 ks XMM-Newton obser- 
vation of NGC 6860. While the spectral state of the source changed between the 
newer observations presented here and the earlier short XMM-Newton spectrum 
- showing a higher flux and steeper power law component - the spectrum of 
NGC 6860 is still complex with clearly detected warm absorption signatures. We 
find that a two component warm ionized absorber is present in the soft spec- 
trum, with column densities of about 10^° and lO^^cm"^, ionization parameters 
of ^ = 180 and 45 ergs cm s~^, and outflow velocities for each component in the 
range of ~ 0-300 kms'^. Additionally, in the hard spectrum we find a broad 

11000 kms~^) Fe Ka emission line, redshifted by ~ 2800 kms~^. 

Subject headings: X-rays: galaxies: active 



1. Introduction 



Through the very hard X-ray survey (14 - 195 keV) conducted by Swift's Burst Alert 
Telescope (BAT), a sample of unbiased t owards absorption (for A^h < lO^^cm"^), local 
(<z>^ 0.03), AGN have been identified ( iTueller et al.l 120081 ). Among these 153 sources, 
we obtained 10 ks XMM Newton follow-ups for 22 previously unstudied, in the X-ray band. 
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BAT AGNs. These sources, selected purely on their significant detection (> 4.8cr) in the 
BAT and probable optical identification with bright sources in the 2MASS/Digital Sky Sur- 
vey, comprise a representative sample of the BAT AGNs, with half of the so urces showing 



low X -ray column densities (iVn < 10^^ cm ^) and half more heavily absorbed (IWinter et al. 



20081 ). From our analysis of both these XMM Newton follow-ups and a combination of 



published spectra and our analyses of S wift XRT/ASCA/XM M Newton spectra of the re 



maining sources in the 9-month catalog (jWinter et al.ll2009al ). we found that the majority 



of the X-ray spectra were well described by simple models, such as an absorbed power law, 
partial covering, soft excess, warm absorber (through fits with O VII and O VIII edges), 
or refiection. Such simple models can even be applied to the more complicated spectra of 
Compton-thick sources, where a strong Fe Ka fiuorescent line and fiat power law spectrum 
indicates a refiection-dominated spectrum. 

However, such simple models could not fit the 0.3-10 keV spectrum of roughly 10% of 
the BAT-detected AGN, including NGC 6860. Unlike any previously known sources, the 
spectrum of this source was fiat (F < 1.0) both below and above 2keV, no matter how 
we modeled the continuum. To our knowledge, there is no other known AGN to exhibit a 
fiat spectrum below 2keV. Likely, this fiat spectrum at low energies is due to complex line 
structure that is unresolved in the ~ 10 ks MOS spectra (the pn spectrum unfortunately 
could not be extracted due to problems with the observation data files (ODFs)). Further, 
we found that the spectrum of this source could be explained by a number of models, in 
particular a warm absorber or double partial covering model fit equally well. This degeneracy 
could not even be broken with the inclusion of the 14-195 keV BAT spectra. 

To unravel the complexities of this unusual source, we obtained a long XMM-Newton 
exposure including both EPIC CCD and RGS spectra. We also include an analysis of the 
archived Suzaku observation of this source, which provides simultaneous coverage of the 
soft, hard, and very hard X-rays through the XIS and HXD PIN detectors. NGC 6860 
is a southern hemisp here {a =20:08:46 .89, 5 =-61:06:00.70) target optically identified as 



a Seyfert 1.5 source (iLipari et al.lll993l ). The source is nearby {z = 0.014884) and hosted 
in a barred spiral galaxy (SBab). Optical emission line diagnostics reveal the source as an 
AGN/s tarburst composite, with the AGN dominating only in the inner 10" in the optical/IR 



bands (iBennert et al.ll2006l ). Thus, the X-ray observations will also likely include contribu- 
tions from both the starburst and AGN. In Section [2] we describe the XMM-Newton and 
Suzaku observations. In Section [3], we include an analysis of the light curves. In Section HJ 
we discuss the spectral analysis. We discuss some of the unusual properties of this source in 
Section O Finally, we summarize the results of our analysis in Section O 
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2. Observations and Data Reduction 

2.1. XMM-Newton Observations 

NGC 6860 was observed on 29 March, 2009 for 123 ks (observation 0552170301). The 
XMM-Newton data were reduced using the Science Analysis System (SAS) version 8.0. 
Calibrated photon event files were created for the EPIC pn and MOS cameras using the 
ODFs and the commands epchain and emchain, respectively. The event tables were filtered 
using the standard criteria outlined in the XMM-Newton ABC Guid^. For the MOS data, 
good events correspond to a pulse height in the range of 0.2-12 keV and event patterns that 
are characterized as 0-12 (single, double, triple, and quadruple pixel events). For the pn 
camera, only patterns of 0-4 (single and double pixel events) are kept, with the energy range 
for the pulse height set between 0.2 and 15 keV. Bad pixels and events too close to the edges 
of the CCD chips were rejected using the stringent selection expression "FLAG = 0" . 

The source and background spectra, along with response and ancillary response files, 
were then extracted using the SAS task especget. The source spectra were extracted from 
a circular region centered on the source, with a radius of 120". The background spectra were 
extracted from a circular region of the same size on the same chip. The spectra were then 
grouped using the FTOOL grppha, binning the spectra by 20 counts. The resultant exposure 
times of the spectra were 100 ks (pn) and 117ks (MOSl and M0S2). The corresponding 
count rates for the spectra are 9.5 counts s~^ (pn) and 2.6 counts s~^ (for MOSl and M0S2). 

Since the XMM-Newton spectrum of our source suffers from high background count rates 
(particularly during the last ^ 8 hours of the exposure), we also extracted a pn spectrum 
from an event file filtered by count rate. We chose the stringent filtering of keeping only 
events whose count rates are below 100. The filtered pn spectrum has an average count rate 
of 9.8 counts s~^ and an exposure time of 66 ks. We compared the spectra of the filtered and 
unfiltered pn count rates. We find that the only differences in spectral shape occur above 
about 8 keV. At these high energies, the unfiltered spectrum shows a curvature which is not 
replicated in the filtered spectrum. The filtered high energy spectrum is instead well fit by 
a power law. Throughout the text, we use the unfiltered spectrum, which has more counts 
and thus a higher signal-to-noise ratio, for all but the broad band spectral fits (for which we 
use the filtered pn spectrum). 

For the RGS spectra, we extracted first-order spectra for the RGSl and RGS2 using 
the SAS command rgsproc. Due to flaring seen in the light curve, we created a good 



http : //heasarc . gsf c . nasa. gov/docs/xmm/ abc/ 



-4- 



time interval file using the SAS task tabgtigen and filtering out all events with count rates 
> 0.2. The spectra were re-extracted using this filtering. Response files were created using 
the task rgsrmf gen and the spectra were rebinned to 20 counts per bin using grppha. The 
exposure times for the RGS spectra are 93 ks, with count rates of 0.21 counts s~^ (RGSl) 
and 0.25 counts s"^ (RGS2). 

2.2. Suzaku Observations 

In addition to the XMM-Newton observation, we include in our analysis the archived 
44 ks Suzaku observation of NGC 6860 (observation 703015010) from 7 April 2008, ob- 
served in HXD nominal pointing mode. We used the cleaned version 2.0 processed event 
files supplied by the Suzaku team. To extract the spectra, we combined the 3x3 and 
5x5 edit modes for each of the front-illuminated (XISO and XIS3) and back-illuminated 
(XISl) CCDs. Source and background spectra were extracted in circular regions of 210"and 
180", respectively, in XSELECT. The response and ancillary response files were created us- 
ing xisrmf gen and xissimarf gen. We then added the spectra and response files from the 
front-illuminated detectors (using mathpha, addrmf , and addarf ) and rebinned the spectra 
and rmf files from 4096 to 1024 channels (which is still larger than the resolution of the 
CCDs). The front-illuminated and back-illuminated spectra were each binned to 20 counts 
per bin. The resultant exposure times are 44 ks with count rates of 1.34 counts s~^ (XISl) 
and 1.13 counts s~^ (the average for the combined XISO -|- XIS3 spectrum). 

For the HXD spectra, we only extracted spectra from the PIN instrument. Spectra 
were extracted using XSELECT for both the source and Suzaku team supplied instrumen- 
tal tuned background file. The source spectrum was corrected for instrument dead time. 
We also included the contribution from the cosmic X-ray background, as suggested by the 
Suzaku teanfl, as a cutoff power law of the form CXB(E) = 9.0 x 10"^ x (E/3keV)-°-29 x 
exp(— ii^/40keV) erg cm~^ s~^ sr^^ keV~^. The standard response file from the Suzaku 
CALDB was used and the spectrum was binned to a signal-to-noise ratio of 3cr. The net 
exposure time for the PIN spectrum is 34 ks with a count rate of 0.31 counts s~^. 



^See the Suzaku ABC Guide, Section 7.5.4 at 
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2.3. The Swift BAT Spectrum 



To fully characterize the spectrum of NGC 6860, we also include an analysis of the Swift 
BAT spectrum generated from the first 22 months of the all-sky survey. The BAT spectrum 
is an 8-channel spectrum in the 14-195 keV band. B oth the spectrum a nd diagonal matrix 
are further described in the 22-month catalog paper (ITueller et al.ll2010f ). 



3. Light Curve Analysis 



In our previous analysis of variability in NGC 6860 ( iWinter et al.l 120081 ). we found 
that the EPIC MOS spectra were variable over the ~ 10 ks observation. We also found 
that a comparison of this short XMM-Newton observation with the Swift XRT observation 
showed a change in the soft (0.5-2 keV) flux with no change in the hard (2-10 keV) flux. The 
difference between these observations was about 3 months. Therefore, the source is known 
to vary in the soft band both on short timescales and over months. 

To determine the nature of variability over the long XMM-Newton observation, we 
constructed light curves for the EPIC pn and MOS spectra of NGC 6860. We constructed 
light curves in three separate bands: soft (0.1-3keV), medium (3-7keV), and hard (7- 
lOkeV). In this way, we can associate the amount of variability with features prominent in 
each of these bands - namely, a soft excess and/or warm absorbers in the soft band, the 
region surrounding the Fe Ka region for the medium band, and the hard continuum in the 
hard band. In Figure [H we plot the background subtracted pn light curves, binned on a short 
1024s timescale and the orbital timescale (5760s, showing the normalized light curves). The 
mos light curves, which are not shown, exhibit the same trends seen in the pn. 

The light curves show that the soft band dominates the total count rate and that this 
band varies by count rates of as much as 2 counts s~^ over the observation. The mean 
count rates in each of the bands correspond to 6.9 counts s~^ (soft), 1.3 counts s~^ (medium), 
and 0.2 counts s~^ (hard). When we look at the normalized light curves, we flnd that the 
variability in the soft band correlates well with the variability in both the medium and hard 
bands. 

In a similar way as for the XMM-Newton light curves, we constructed light curves in a 
soft, medium, and hard band for the Suzaku XIS detectors. Since the response for the XIS 
detectors extends to higher energies, we used the energy range of 7-12 keV for our XIS hard 
band. We created a combined (XISO + XISl + XIS3) background subtracted light curve for 
the Suzaku observation, binned by 128 s. As for the XMM-Newton observation, we found that 
the low energy band had the highest average count rate (< Ct Rate >= 0.83ctss~^). Back- 
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ground subtracted mean count rates in the medium and hard band correspond to 0.34ctss~^ 
and 0.05 ctss^^, respectively. In Figured we plot the normalized 128 ks light curves. As for 
the XMM-Newton observations, we find that the variability in the soft, medium, and hard 
bands is well-correlated, particularly for the soft and medium bands. 



4. Spectral Analysis 

Given the complex nature of the spectra of NGC 6860, we began our spectral fits by 
analyzing the spectra in the soft and hard energy bands separately. We first determine the 
parameters of the hard X-ray continuum and the region surrounding the Fe Ka emission 
in Section 14. 1[ This analysis includes spectra from Suzaku and XMM-Newton, focusing on 
the spectral region from 3-50 keV. To determine the properties of the soft band, the region 
corresponding to the highest count rates for NGC 6860 and the greatest complexity, we 
examine the XMM-Newton and Suzaku CCD spectra as well as the XMM-Newton RGS 
spectra in Section 14.21 We then fit the broadband continuum, including the Swift BAT 
spectra, in Section 14.31 We use the best broad band fit to constrain the continuum in the 
RGS data, to better determine the properties of the warm absor bers in Section 14.41 All of the 



spectral fits described were done using XSPEC version 12.5.1 (lArnaudlll996l ). The quoted 
errors correspond to a change in = 2.71 or the 90% confidence level for an additional free 
parameter. 



4.1. The Hard X-ray Continuum and the Fe Ka Emission Region 

In the following sub-sections, we describe spectral fits to the hard X-ray emission 
(> 3keV) for the Suzaku and XMM-Newton observations. We include absorbed power law 
models, as well as fits with a cutoff energy/reflection component. Also, we fit for additional 
lines/absorption edges in the region around the 6.4 keV Fe I Ka emission. 



4.1.1. The Suzaku spectra 

We conducted combined fits of the Suzaku XIS (XISl and XIS0+XIS3) and PIN spectra, 
constraining the PIN data such that its spectrum is fixed at a constant level of 1.16 of the 
XISl spectrum, as indicated in the Suzaku Data Reduction Guide§|. We find that a simple 
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cm ^ with a multiplicative tbabs model (IWilms et al.l |2000| ) , we find that the neutral 



cm ^ with the multiplicative 



absorbed power law model provides a good fit to the dat a (red uced x^/dof = 1.09/1085) 
Fixing the Galactic absorption to the Dickev fc LockmanI (199o[) measured value of 4.19 x 
10 

hydrogen colum n density of the sou rce is best fit as 4.8 

ztbabs model (jWilms et al.l 120001). The power l aw component measured is typical of the 
Swift BAT-detectected AGN Jwinter et al.ll2009ah . with T = 1.70±0.05 (from the zpowerlw 
model). Residuals to the model indicate the presence of an Fe Ka emission line. 

Adding a Gaussian to account for the Fe Ka emission with the zgauss model, we 
find that the line is significant with an improvement in the fit by Ax^ = —48. The mea- 



sured parameters for this line are E — 6.31 ± 0.05 keV, a — 0.15^o;o4 keV, / 



2.6 



+0.6 



10 ^photons cm ^s ^, and EW = 102l32eV. The width of this line indicates an origin in- 
terior to the broad line region (with AGN H/3 FWHM values typic ally ranging from 30 00 
- 10,000 kms"^ and measured as 5920 ± 600 kms"^ for NGC 6860 feennert et aDl2006h ). 
with a corresponding velocity of ~ 17000 km s~^. This broad width is statistically a better 
fit than a narrow line fit, with a Ax^ = —7.6 over a fixed width of cr = 0.01 keV. 

The addition of a 7.11 keV Fe edge (over the edge built into the ztbabs m odel), is also 
signif icant, with Ax^ = —6.1. The energy of the edge is consistent with cold Fe ( iHenke et al. 
19931 ). with E = 7.11^Q;j^g keV. The measured optical depth corresponds to r = 0.08 ± 0.05 



The Fe edge can be the result of either reflection of the direct emission, a higher than solar 
iron abundance, or a higher column density of gas than measured fro m the neutral absorp tion 
model. Computing the Fe XVII cross section at 7.11 keV following IVerner et al.l (Il996[ ). we 
find a cross section of cr{E) = 4.5 x 10~^Mb. From the estimated cross section and our 
measured optical depth, we find a column density of iron that is 1.80 x lO^^cm"^ (iVpe = 
T/a{E)). Therefore, the estimated iron abundance is very high (A^pe/A'^H ~ 139(Fe/H)0) if 
the edge is associated with the 4.8 x 10^^ cm~^ absorber measured in our fits. If, however, 
the abundance of iron in NGC 6860 is the solar value, this implies that the Fe edge is 
instead the result of a much higher hydrogen column density of ~ 7 x 10^^ cm~^, which 
does not affect the emission in the softer bands (i.e. a partially covering near Compton- 
thick absorber). Finally, while the residuals to the model suggest there may be additional 
emission/absorption features present, we found no statistical support for additional emission 
features found in the spectra of many AGNs (Fe XXV Ka, Fe I K/3, Ni I Ka). 

In addition to a simple absorbed power law model, we tried to determine whether a high 
energy cutoff or reflection component are present in the data. To constrain a possible high 
energy cutoff, we substituted the zpowerlw model with a cutoff power law model. However, 
we found no measurable improvement in the fit and could not constrain the energy of the 
cutoff. 
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Additionally, we tested for the presence of a reflection component, which could alter- 
natively account for the presence of the 7.11 keV Fe edge. Using the pexrav model for 
reflection off of cold neutral material, we tied the normalization and power law index to 
be the same values between the direct power law component and the reflection component. 
Statistically, the reflection model does not significantly change the fit, improving by only 
4.4, while losing two degrees of freedom. Using this model, we find neither the folding energy 
{Efoid = 30l^9 keV) nor the reflection parameter (which can range from 0-1 with no change 
in x^) are well constrained. These fits, however, do not include the higher energy BAT data, 
which are used to place tighter constraints on reflection and the energy cutoff in Section 14.31 

4.1.2. The XMM- Newton EPIC spectra 

For spectral fits of the hard continuum in the EPIC spectra, we simultaneously fit the pn, 
mosl, and mos2 spectra in the 3-10 keV band, using a constant value to allow for variations 
in the flux levels between CCDs. We find a good fit (reduced x^/dof = 1.17/2336) with a 
simple absorbed power law model. The parameters of the model indicate a column density 
(above the Galactic value) and power law index similar to those measured in the Suzaku 
observation: A^h = ^.Slil'l ^ ^^'^^ cm~^ and F = 1.69 ± 0.04. The residuals to this model, 
shown in Figure [31 clearly show an Fe Ka emission line. Additionally, residuals suggest an 
emission feature red- ward of the neutral Fe Ka emission and possibly an 7.11 keV Fe edge (in 
the pn spectrum). Fitting an Fe edge at 7.11 keV, however, results in no significant change 
in of the model. We find an upper limit on the optical depth of the edge at r = 0.07, 
which is in agreement with the Suzaku results. 

The addition of a Fe Ka emission line is clearly significant, resulting in Ax^ = —190.6. 
The measured line parameters are again similar to those found in the Suzaku observation. 
We find E = 6.35 ± 0.02 keV, a = Q.Uto^ol^eV, I = 2.21°;^ x 10"^ photons cm'^ s-\ and 
EW = 85li\leV. The line is clearly broad, since there is a Ax^ = 15 between the fitted 
broad a and a fixed narrow line (a = 0.01 keV). The best-fit energy, width, and equivalent 
width are consistent between the XMM-Newton and Suzaku observations. Since we would 
expect an Fe K/3 emission line to also be present if the Fe Ka emission is produced in 
neutral material, we added this feature to our spectral fits. We fixed the centroid energy to 
7.056 keV, the width to 0.01 keV, and the intensity to be 13% of the Fe Ka line. There is 
virtually no change in x^ (Ax^ = —0.01) and no change in the parameters of this fit upon 
adding the Fe K/3 line. 

Adding a narrow Gaussian (cr = 0.01 keV) to fit the additional emission line seen in 
the residuals, we find that the line is significant, improving x^ by 15. The fitted parameters 
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of this emission feature are E 



5.951°:°^ ^eV, / 



5.1^^3 X 10 ^photons cm ^s ^, and 



EW = 17lg eV. While the hne is weak, it is also significant. The energy of the line is unusual, 
but c orresponds to line s seeii in the spectra of other AGN sources. For instance, as pointed 
out in lMarkowitz et al.l (120091 ). high resolution spectroscopy of 16 Seyferts have shown narrow 



features interp reted as blueshifted or r edshifted Fe K features from hot spots above the 
accretion disk (IVaughan fc Uttleyll2008l ). The hne we have detected has roughly the same 
energy , width, and strength as the 6.0 keV feature detected in NGC 3227 by lMarkowitz et al. 
J2009h . 



Comparing the flux from the best fit absorbed power law + Fe Ka emission model 
with XMM-Newton to the Suzaku spectra, we find no change in flux. The XMM-Newton 
observation has a 3-10 keV observed flux of 1.86 x 10^^^ ergss^^ cm~^ in the pn. The XIS 
spectra have an average flux of 1.84 x 10^^^ ergs s^^ cm~^. Therefore, while the XMM-Newton 
observations were taken a year after the Suzaku observations, there is virtually no change in 
the hard X-ray spectra. 



4.2. The Soft X-ray Spectra 

In the following sub-sections, we detail simple model fits to the soft spectra from the 
Suzaku and XMM-Newton observations. The base model is an absorbed power law, at the 
redshift of the source. To account for the additional features in the soft spectra, we use a 
combination of absorption edge models and Gaussians to fit absorption and emission features. 



4-2.1. The Suzaku XIS spectra 

Fitting the soft XIS spectra simultaneously in the 0.3-3 keV band, we find that a simple 
absorbed power law model is not a good fit to the data with a reduced x^/dof = 1.64/883. 
Comparing the fitted parameters of this model to the hard continuum, we find that both 
the neutral hydrogen column density (A^h = 2.9 ± 0.9 x 10^°cm~^) and power law index 
(r = 1.36 ± 0.03) are lower in the soft band. As shown in Figure SJ the simple power law 
model clearly can not account for the soft spectrum. 

We find that the addition of warm absorption signatures, through O VII and O VIII 
edges, to the model greatly improves the fit. Figure S] shows the improvement in the fit 
for the ratio of data/model. With these edges added to the model, we find Ax^ = —514.7 
and a reduced x^/dof = 1.07/879. Additionally, while the column density of the intrinsic 
absorption is still low (A^h = 7.9to9 x lO^'^cm"^), the power law index from this warm 
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absorber model is consistent with the hard power law continuum (F = 1.74 ± 0.04). The 
best-fit parameters for the oxygen edges correspond to E = 0.73 ± 0.01, r = 0.34 ± 0.04 
(O VII) and E = 0.87 ± 0.01, r = 0.35t[}:|j^ (O VIII). In addition to the O VII and O VIII 
edges, the residuals to the model suggest another edge feature near 2keV. Adding an edge 
improves the fit by Ax^ = 10. This edge, which has best-fit parameters of £^ = l.Ol^l^gogkeV 
and r = O.Obbtomr, 

is likely associated with errors in the detector calibration of the Si 
K edge. For our final best-fit model, the estimated observed flux in the 0.3-3 keV band is 
1.34 X 10^^^ ergss^^ cm"^. 



4.2.2. The XMM- Newton EPIC spectra 

Given the much higher signal-to-noise in the XMM-Newton EPIC spectra, a simple 
absorbed power law model does not adequately describe the 0.3-3 keV EPIC spectra, with 
a reduced x^/dof = 7.71/898. In Figure El we plot the ratio of the data/model for the 
simple absorbed power law model. As for the Suzaku spectrum, the best fit requires a lower 
neutral hydrogen column density (A^h = 2.4 x 10^°cm~^) and a flatter power law continuum 
(F = 1.46). However, as shown in the figure, there are more features than can be accounted 
for with the addition of O VII and O VIII edges. 

While adding the O VII and O VIII edges improves the fit significantly, with Ax^ = 
—4458.0, the fit is still not acceptable (reduced x^/dof = 2.75/894, see Figure[5]). With this 
higher signal-to-noise data, a better interpretation for the residuals between 0.7-1.0 keV, from 
the absorbed power law model, is that this feature is from a broad unresolved transition array 
(UTA) from Fe M-shell ions. Such features are detected in the XMM-Newton and Chandra 
spectra of a number of Seyfert 1 sources (e.g. NGC 3227, NGC 3783, NGC 4051, Mrk 509). 

We find that the addition of Gaussians/edges to the data greatly improves the fit. 
However, we note that due to inconsistencies in the pn and MOS calibration below 0.5 keV, 
we used only the pn data to determine properties of lines detected in this low energy region. 
In Table [T], we specify the parameters of lines/edges where |Ax^| > 15. These features 
include N VI, Ne IX, and Si XIII absorption lines. We found the addition of another edge 
model, at 1.7keV, improved the fit greatly (Ax^ = —26). The most probable identification 
for this feature is the Mg XI edge with a rest energy of 1.761 keV, which is not near any 
known energies corresponding to calibration issues for the EPIC pn. The improved fit results 
in a reduced x^/dof = 1.45/855. 

From 0.45 to 0.5 keV, the residuals showed complex absorption/emission that could not 
be resolved in the pn. We found that an emission line at 0.45 keV greatly improved the 
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fit. Given the uncertainties in the pn cahbration below 0.5 keV (see the EPIC cahbration 
document0), we can not uniquely identify this feature at 0.45 keV. However, we note that 
this is possibly N VII. 

The associated outfiow velocities determined from the fitted energies of the edges and 
lines are difficult to interpret. We find that the velocities of the O VII and O VIII edges are 
consistent with each other, implying red-shifted velocities of approximately 8000-9000 km s~^. 
The estimated velocity of the Mg XI edge is consistent with the oxygen edges. However, the 
energy of the Ne IX emission line suggests a blue-shifted velocity of about 5000 km s~^ while 
the remaining lines, are consistent with no velocity shift. The inconsistent velocities are 
likely due to a combination of poor energy resolution in the EPIC spectra and confusion of 
multiple features. The 20-30 eV shifts in the fits of the spectral features seen in the pn data 
are probably due to instrumental effects, since they are not seen in the Suzaku data. Shifts 
this large, while rare, have been seen in pn, particularly during episodes of high background. 

The best fit neutral column density is much lower than that found from fits to the 
hard spectrum (we find Nh = 5.93lo'g7 x 10^°cm~^). However, we find that the best-fit soft 
band power law index is consistent with the hard band value (we find F = 1.697^qqq3). The 
estimated fiux in the 0.3-3 keV band is 1.31 x 10^^^ ergss"^ cm~^, which is also consistent 
with the soft fiux in the Suzaku observation. 



4.2.3. The XMM-Newton RGS spectra 

As a first fit to the RGS spectra, we fit an absorbed power law model (zpow), with 
absorption fixed at the Galactic value. Unlike with the EPIC CCD data, we find a good fit 
(x^/dof = 2865.9/1960) with a slope, P = 1.50 ± 0.02. The addition of an intrinsic neutral 
absorber does not improve the fit. The addition of a blackbody component improves the fit 
significantly, with Ax^ = —105.3. Adding edge models to account for the O VII and O VIII 
absorption edges further improves the fit with Ax^ = —129.3. The best-fit parameters 
are: O VII edge energy = 0.705toml^eV with optical depth r = 0.15+[^:[!^, O VIII edge 
energy = 0.850lJ^:[!|J^ keV with optical depth r = O.SOtoM, and kT = 0.30lJ^:[!^ keV with 
normalization of 1.9lg7 ^ 10^^ photons cm'^ s^^ at 1 keV for the blackbody component. The 
fiux in the 0.3-2.5keV band is 1.22 x 10"^^ ergscm"^ s~^. 

Following this simple power law + blackbody fit (along with the O VII and O VIII 
edges), we determined the significance of additional absorption and emission features seen 
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in the residuals to the model. We significantly detect (i.e. |Ax^| > 10 when a Gaussian 
is added) N VI, O VII, and Ne IX in absorption and N VII in emission. Details of these 
Gaussian fits, including the energy and intensity of the lines, are shown in Table [2J For all 
of these fits, we fixed the width of the lines to 0.01 keV. An additional absorption line of 
Mg XI appears in the residuals, but is not as significant in A^^. We show the ratio of the 
data/model for the RGS spectra in Figure [61 In addition to the features mentioned, there 
is evidence of an Fe M-shell UTA between 0.7 and 0.8 keV. While residuals to the model 
suggest that a Mg XI edge may be present in the RGS spectrum, adding the feature does 
not statistically improve the fit. 

The combination of lower signal-to-noise in the RGS spectra and worse energy resolu- 
tion in the EPIC spectra make it difficult to determine why some lines are only detected 
significantly in either the RGS or EPIC spectra. However, with the exception of the Si XIII 
line, the discrepancy in detection of lines is at the low energy range of the bandpass, from 
about 400-550 eV. In this region, there are uncertainties in the pn calibration. Additionally, 
there are many features which are not easily resolved even with the higher resolution in the 
RGS spectra. Therefore, we can not fully determine the emission and absorption properties 
of lines in this region without a higher signal-to-noise spectrum with the energy resolution 
of a grating spectrometer. The lack of a detection of Si XIII absorption in the RGS could 
either be the result of the lower signal-to-noise of the RGS (the intensity of the line may be 
too low to detect significantly in the RGS spectra) or this line may be due to calibration 
uncertainties associated with the Si K edge in the EPIC detectors. One clear conclusion 
we can draw from the comparison of the EPIC and RGS data is that the most significant 
features are the O VII and O VIII edges and the Ne IX absorption line. 

Unlike our analysis of the soft EPIC spectra, we find that the implied velocities from the 
oxygen edge energies do not agree with each other. Assuming the lines are resonant lines, 
the energies of the N VI and O VII lines are consistent with each other, as well as the O VIII 
edge, with blue-shifted velocities of about 7000 km s~^, similar to the velocities determined 
for the edges measured in the EPIC data. The energy of the O VII edge suggests a velocity 
twice this value. However, the energy of this edge is difficult to determine due to confusion 
with the Fe M-shell UTA (as seen in Figure E]). Meanwhile, the energies of the Ne IX and 
N VII lines suggest low to no velocity shift. These results suggest that there may be two 
different velocity components in the ionized absorber. We further investigate the properties 
of the ionized absorbers with warm absorption model fits in Sections 14.2.41 and 14.41 
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4 -2. 4- Warm Absorber Model Fits 

Having fit the spectra of our sources witli simple models for tlie absorption/emission 
lines in the soft spectra, we next investigate the properties of the warm absorbers through use 
of an analytic model. The model warmabs is an analytic version of the XSTAR table models 



(iKallman &: McCraylll982l ) available as additional models to XSPEC. The model calculates 
spectra using stored level populations, which are scaled using the elemental abundances 
specified by the user. The free parameters in the model include the column density of the 
absorber, the ionization parameter (log^, where ^ = L/(ni?^) and ^ has units of ergss"^), the 
redshift of the ionized gas, and the turbulent velocity of the gas (which affects the broadening 
of the spectral features). For the warmabs model, we assume a gas density of 10^ cm~^ and an 
illuminating power law with F = 2, as indicated in the population model pops . fits.n04- We 



assume solar abundance for the elemental abundances, using the abundances of IWilms et al. 



(I2OOOI ) (in XSPEC, abund wilm). Additionally, we assume an outflow velocity of kms~ 
for these initial fits, fixing the redshift of the warm absorber to the source redshift. In 
Section 14.41 we better constrain the continuum in the RGS spectra and re-fit the ionized 
absorbers allowing the velocity of the outflow to vary. 

We used this analytic warm absorption model to determine the best-fit parameters to 
the soft X-ray spectrum (0.3-3 keV). First, we fit each of the spectra individually (Suzaku 
XIS03+XIS1, XMM pn, and XMM RGS) with a base power law model, absorbed by the 
fixed Galactic absorption. Since the RGS data required a black body fit, an indicator of 
the 'soft excess', we added this as well as a component for neutral intrinsic absorption to 
the model fits. The residuals to all of the spectra showed a poor fit, particularly to the 
Fe UTA feature. Next, we added a single ionized absorber with the warm absorber model. 
In all cases, the addition of a warm absorber model significantly improved the model fit 
and provided a good fit to the Fe UTA feature. We then added a second warm absorber 
model, to determine if a second component is required. For both the Suzaku XIS and XMM 
RGS spectra, the addition of a second component is statistically significant (|Ax^| > 50, see 
Figure El for the ratio of the data/model for the RGS spectra). For the XMM pn spectrum, 
a second component is not required. However, the residuals to the model for the warm 
absorption models show a poor fit below ^ 0.7 keV. The best fit parameters to the soft 
spectra with the warm absorption models are shown in Table O 

From our neutral + warm absorber model fits, we find that the neutral column does not 
change drastically between the Suzaku and XMM-Newton observations. Additionally, the 
RGS and pn columns are fairly consistent (Nh ~ lO^^cm"^). Comparing the XIS and pn 
observations, we find that the warm absorber (component 1 in the XIS) is consistent within 
the measured errors. These values point to an ionized absorber with a column density of 
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^ lO^^cm"^, an ionization parameter of log.^ ~ 2.0, and a turbulent velocity from ^ 25- 
200kms~^. The most notable differences between the Suzaku and pn spectra include the 
normalization on the black body component (which is roughly twice as high in the XIS 
spectra) and the slope of the power law (which is flatter in the pn observation). 

There is clearly a mis-match between the best-fit parameters of the XMM pn and RGS 
spectra, which were taken simultaneously. The inconsistencies are most apparent in both 
the blackbody (the RGS value is ~ IkeV, while the pn value is typical of the soft excess 
with kT = 0.13keV) and power law component (the RGS value is too steep, while the pn 
value is too fiat). With the warm absorber model, both the reduced (1.25 for RGS and 
1.36 for pn) and residuals to the model show that the fits are not perfect representations of 
the spectra. Further, the higher signal-to-noise in the pn and the better energy resolution 
in the RGS make it difficult to reconcile discrepancies between the fits. 



4.3. Broad Band Spectral Fits 

Since we found no significant difference in either flux or spectral shape between the 
Suzaku and XMM-Newton observations, we fit the Suzaku XIS and XMM-Newton pn data 
simultaneously with the Suzaku PIN and Swift BAT spectra to obtain the best broad-band 
(0.3-195 keV) spectral fit. A major problem of fitting non-simultaneous data, like the XMM- 
Newton pn and Swift BAT spectra, is constraining the proper normalization between the 



softer and hard energy spectra (see IWinter et al.l (120081 ) for a discussion of problems with 



joint XMM-Newton and Swift spectral fits). This problem is mitigated by the addition of the 
Suzaku data, where the normalization between the PIN and XIS is known. A second potential 
problem with using the Swift BAT spectra is that the BAT spectra are time averaged over 
a 22-month period. Therefore, any changes in the spectral shape and/or flux will not be 
accounted for in the time averaged spectrum. For NGC 6860, however, the Swift BAT 
54- month light curve shows no significant variation (private communication, Gerry Skinner). 

For the broadband spectral fits, we used a constant value (const) to account for dif- 
ferences between the Suzaku XIS and the additional spectra. This constant was set to 1.16 
for the Suzaku PIN spectrum but allowed to vary for the XMM-Newton pn and Swift BAT 
spectra. Our base model is a power law affected by the Galactic neutral column density. To 
account for the warm absorber features, we added edges to account for the O VII and O VIII 
absorption edges and Gaussians to account for the N VII and Ne IX absorption lines found 
in the XMM-Newton pn and RGS spectra. A blackbody model was added to account for 
the soft excess. Also, we added a Gaussian to account for the Fe Ka emission near 6.4 keV. 
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The addition of an intrinsic absorption model (with the ztbabs model) to the spectral 
fits greatly improves the fit with Ax^ = —3025.8. The fitted column density, Nh ~ 6 x 
lO^^cm"^, is consistent with the neutral column derived from the warm absorber model fits 
to the soft spectra. To determine if the absorber is partially or fully covering the spectrum, 
w e tried substit u ting th e fully covering ztbabs model with the pcf abs model. Previously, 



m 



Winter et al.l ( l2009al ). we found that the partial covering model was successfully applied 
to the majority of complex spectral sources in the Swift BAT AGN catalog. The partial 
covering absorption model (pcf abs) has two parameters, the covering fraction (/) and the 
neutral hydrogen column density (Nh), such that: M{E) = / exp(— NHcr(i?)) + (1 — /). 
The parameter o"(E) is the photoelectric cross section. Upon substituting in this model, we 
find that the covering fraction is not well constrained, with a range of possible values from 
0.0-0.96. Additionally, there is no change in when the pcf abs model is used in place of 
the ztbabs model. 

Adding an energy cutoff model to the power law component improves the fit significantly 
(Ax^ = —25). In this model (cutoffpl), the power law component is multiplied by a 
factor of M{E) = exp{—E/Ec), where Ec is the cutoff energy. The best fit cutoff energy 
is lOSlg^ keV, which is consistent with a Comptonization spectrum. Best fit parameters for 
this full model (in XSPEC, the model is ztbabs*zedge*zedge*((3 zgauss components) + 
bbody + cutoff pl)*const) are shown in Table H] (we do not include the parameters for the 
soft energy edges and lines, instead see Tabled]). Additionally, we plot the best-fit spectrum 
in Figure [Si 

To determine the importance of reflection in the spectrum of NGC 6860, we i ncluded re- 



flecti on with the XSPEC pexrav model for Compton reflection off neutral material (IMagdziarz fc Zdziarski 



19951 ). We fixed the normalization and power law component to be the same between the 
reflection and cutoff power law models. We also fixed the folding energy of the reflection 
component to be the same as the power law cutoff energy. In the pexrav model, we con- 
strain the inclination of the disk to be 45° and the Fe abundance to be the solar value. The 
best-fit reflection model is statistically worse than the cutoff power law model (Ax^ = 40), 
but still an acceptable fit to the data with a reduced value of 1.11. We find that the 
best-fit reflection parameter [R = Q/2n, where fl is the angle subtended by the reflector) 
is R = 1.15^Q 4g. With this reflection model, the best-fit neutral column density is slightly 
lower {Nn = G.lStoH x lO^^cm^^). The power law index (F = 1.66+0:04) and cutoff energy 
{Ecutoff = 56+^5 keV) are slightly different than in the direct power law model. The relative 
normalizations from this model (relative to the Suzaku XIS spectra) are the pn at 1.09, the 
BAT at 0.94, and the PIN again fixed at 1.16. 
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4.4. Better Constraining the Properties of the Warm Absorber 

Given the small band pass for the RGS spectra, the AGN continuum emission can not 
be constrained without the use of higher energy data. In our earlier fits of the RGS spectra, 
we found that the fitted power law continuum was much steeper than the pn in a similar 
soft bandpass. The value from the soft band pn continuum fits was also not well-constrained 
- with a spectral index much flatter than that of the broad-band continuum fits. Thus, to 
develop the best constraints on the warm absorption, we used the best-fit power law value 
from the broad-band fits along with the high resolution RGS spectra. We found it necessary 
to fix these power law parameters in the model fits, due to degeneracies between the power 
law, soft excess, and absorption models in this narrow bandpass. 



In the spectral analysis of a similar warm absorber source, NGC 3227, iMarkowitz et al. 



(120091 ) compared three different models to the soft excess in a broad band EPIC spectral 
fit: a steep power law, a blackbody, and a Comptonized component (using the CompST 
model). The best statistical fit was with the Comptonization model (though it was only 
slightly better than a steep power law model). Therefore, we chose to use this model in 
our new RGS spectral fit. The values of the parameters from our best-fit model for the 
Comptonization parameter are a temperature of kT = 1.09 keV, an optical depth of r = 12, 
and a normalization of 1.7 x 10^^ phkeV"^ cm~^ s~^. We could not constrain the errors on 
this model as the parameters of this Comptonization model were not well-constrained. 

The final model for this RGS spectral fit includes neutral absorption and two warm 
absorber components. To constrain the velocity of the outflow, we allowed the redshift of 
the warm absorber components to vary. However, we kept the turbulent velocity fixed at 
the default value of lOOkms"^. We include the best-fit absorption properties for this model 
in Table |5l Both the measured column densities of the cold neutral and warm absorbers and 
the ionization parameters are well matched between the previous model fits in Table |3] and 
the new values in Table |5l However, we have placed tighter constraints on the model in our 
new fits by using a fixed power law component as determined from the broad-band fits. Also, 
we have placed estimates on the velocity of the outfiow. The associated systemic velocities 
of both components (calculated from the redshifts recorded in Table |5]) are vi = — 184(1^Jq) 
kms~^ and t>2 = — 196(1}qo)- This places constraints on the outfiowing gas to be from about 
0-300 kms^^ 
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5. The Unusual Properties of NGC 6860 



In IWinter et al.l (120081 ) . we pointed out the unusual spectral properties of NGC 6860 
based on a ~ 10 ks XMM- Newton follow-up observation. The flat (F ^ 1.0) spectrum at hard 
energies led us to classify the source as a Compton thick candidate. Compton thick sources 
are those where the optical depth towards Compton scattering is r >> 1, corresponding to 
column densities A^h ^ lO^'^cm"^. Compton thick spectra are reflection dominated, since 
little to no direct emission escapes below 10 keV. Thus, both a flat spectrum and high EW Fe 
Ka emission line - signatures of a reflection spectrum - are diagnostics used in classiflcation 
of a source as Compton thick. Another possible diagnostic is a lack of variability, as discussed 



m 



Winter et all (l2009bf ). 



We classifled the short XMM-Newton observation as a candidate Compton-thick source 
due to the unusual properties of the spectrum. While the source showed a flat spectrum at 
high energies, the measured column density indicated a column density of only ^ 10^^ cm~^. 
Additionally, the flat spectral index measured at low energies {E = 0.3-2keV), suggested 
that the source spectrum also includes a contribution from warm ionized gas. Without a 
higher signal-to-noise observation we could not further determine the X-ray spectral prop- 
erties of NGC 6860. 

From our analysis of the new higher signal-to-noise spectra presented in this paper, 
we flnd that the spectra of NGC 6860 do not meet the criteria typically used to classify 
a source as Compton thick. The continuum at high energies is no longer flat - marking a 
clear change in spectral state from the earlier XMM-Newton spectrum. Additionally, the Fe 
Ka line has an equivalent width too low {EW ~ 100 eV) to suggest a reflection dominated 
spectrum. Further, the measured column density is well below the Compton thick threshold. 
Therefore, the source is most likely in a Compton thin state throughout both the Suzaku and 
long XMM-Newton observations. Still, it is unusual to flnd a Compton thin source whose 
flux and spectral shape remain constant over a period of a year, as we have found from our 
comparison of the Suzaku and XMM-Newton spectra. Further yet, the Swift BAT 54-month 
light curve shows no variability at the highest energies over the entire 4.5 year period (private 
communication, Gerry Skinner). 

This lack of variability is surprising considering that the previous analysis of the short 
XMM-Newton and Swift XRT spectra showed variability in flux and spectral shape between 
the two observations spaced only 4 months apart. In order to compare the spectra from the 
high (long exposure XMM-Newton) and low (short exposure XMM-Newton and Swift XRT) 
signal-to-noise observations, we flt the 0.3-10 ke V band of the new E PIC spectra with the 



partial covering absorbed power law model used in I Winter et al.l (120081 ). This flt is not a good 



representation of the high quality EPIC spectra (with a reduced value of 3.4), but it allows 
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us to compare the fitted parameters and fiux between the old XMM, XRT, and new XMM 
observations. The best-fit parameters to the new data are: Nh = 2.7 x lO^^cm"^, with a 
covering fraction of 43%, and F = 1.58. As stated earher, the spectral shape is much different 
than the old XMM observation from 3 years earlier (Nh ~ 4.5 x lO^^cm"^, with a covering 
fraction of < 70%, and F 1-0), but similar to the values from the XRT observation. 
Unfortunately, without a higher quality spectrum during the spectral state exhibited in 
the earlier XMM-Newton observation, we can not confirm that the source changed from a 
Compton thick to Compton thin phase. 

Another noticeable change between the short XMM-Newton and XRT observations and 
the later Suzaku and long XMM-Newton observation is the change in fiux. The Suzaku 
and long XMM-Newton observations are 3 times brighter than the short XMM-Newton 
observation and 60% brighter than the Swift XRT observation. It is unclear why the more 
recent observations do not change - whether we happened to catch the source in a similar 
state or the source somehow changed such that the X-ray spectrum is now constant. However, 
given the history of past variability as well as the low Fe Ka EW, it is unlikely that this 
object is truly Compton thick. Although, it is possible that the source changed from a 



Compton thick state - as is the case for the 'changing- look' AGN (IMatt et al.l 120031 ). 



When we examined the properties of the warm absorbers in NGC 6860, we find more 
evidence that its spectrum is unusual. We find that there are likely absorption and emission 
components in the spectrum, but the only strong, clearly identifiable features at soft energies 
are the Fe M-shell UTA and the Ne IX absorption line. We are unaware of any other 
warm absorber spectra that show Ne IX as the most prominent line in the soft spectrum. 
While other absorption and emission features are present, it is impossible to clearly separate 
the features which are likely blended at the resolution available in the long XMM-Newton 
spectrum. 

One of the clearest indications of this complexity is the discrepancy we get between the 
estimated velocities of the edges and lines identified in the FN and RGS spectra. We find, for 
instance, that the fitted energy of the N VII absorption line shows no velocity shift, while the 
N VI and O VII lines are shifted by an apparent blue-shifted velocity of about 7000 km s~^. 
To add to the complexity, we find that the O VII and O VIII edges may either not be shifted 
or blue-shifted depending on the continuum and fitted energies of nearby lines. Even in the 
RGS, we appear to have insufficient resolution to clearly identify all of the features. Further, 
the warm absorption models, while they give an adequate fit, clearly do not fully account 
for all of the features we find (for instance, the Ne IX line has a higher EW than that fitted 
by the warm absorption models). In order to get a better understanding, particularly of the 
kinematics of the gas, it will be imperative to obtain a UV spectrum of NGC 6860. 
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Finally, another unusual property of the NGC 6860 spectra involves the fluorescent Fe 
line in the hard spectrum. The average measured energy for the 6. 4 keV Fe Ka line from 
Chan dra grating observations of a sample of AGNs is 6.399±0.003 keV (lYaqoob fc Padmanabhan 



20041 ) ■ We, however, measure a central energy of 6.331]'^q;q225 ot red-shifted by about 2800 km s ^ 



To conflrm the measured energy of this line, we replaced the broad Gaussian model initially 
used to flt Fe Ka in the broad-band spectral flt (Table H]) with a more physical model of li ne 
emission from a relativistic accretion disk (using the diskline model of lFabian et al.l (119891 )). 
In this model, we used the default emissivity index (/3 = —2) and outer radius (1000 GM/c^). 
We allowed the line energy, disk inner radius, and inclination to vary. The best-flt diskline 
model, with x^/dof = 4282.4/3872, is statistically a worse fit (A^^ ^ 28) than the simple 
Gaussian. However, there is little noticeable difference in the ratio of the data/model be- 
tween the Gaussian and diskline models, as shown in Figure O The best-fit from the diskline 
model confirms the energy measured from the Gaussian model, with E = 6.327]'^qq}7 keV. 
Additional fitted parameters for the diskline model are an inner accretion disk radius of < 23 
GM/c^, inclination of < 15°, and a normalization of 2.4lQg x 10~^ photons cm~^ s~^ (with 
a line equivalent width of 98 eV). With the diskline model, however, the 5.95 keV line is no 
longer significant in the spectrum. 



6. Summary 

Following upon the analysis of a short ^ lOks XMM-Newton observation of the Swift 
BAT-detected Sy 1.5 NGC 6860, we obtained and analyzed higher signal-to- noise observa- 
tions of this source with Suzaku (~ 40 ks) and XMM-Newton 100 ks). Though the new 
observations show NGC 6860 in a different spectral state - with a higher flux and steeper 
power law index - these higher quality spectra conflrm the complex nature of the X-ray 
spectrum of NGC 6860. Even with the high signal-to-noise in these observations, we flnd 
that the X-ray spectra are still too complex to completely characterize. 

Still, a number of conclusions can be drawn from our analysis. Among these, we flnd 
that while short term variability is observed in both the Suzaku and XMM-Newton CCD 
spectra, both the flux and spectral shape remain similar between these observations though 
they were taken a year apart. The overall spectral shape of NGC 6860 shows clear signs of a 
warm ionized absorber at soft energies along with a soft excess. While a combination of two 
warm absorber components and a neutral absorber along with a base blackbody -|- power 
law model provides an acceptable spectral flt in the lower signal-to-noise of the Suzaku XIS 
spectra, the higher signal-to- noise/resolution XMM-Newton pn and RGS spectra are more 
difficult to describe. Particularly, while the Fe UTA features are well flt with the warm 
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absorber models, we find additional absorption features (like the N VII absorption line) 
which are not well- fit by this warm absorber model below ~ O.TkeV. 

The properties of the warm absorbers in NGC 6860 are consistent with those of other Sy 



Is, as determined thr ough analyses of gratin g spectra (IBlustin et al.ll2005l : iMcKernan et al. 



20071 ). For instance, iMcKernan et al.l ( 120071 ) found in a sample of 15 AGN with Chandra 
HETGS spectra, that type 1 AGN typically have multiple warm absorber components with 
ionizing columns from 10^°-10^^ cm~^, ionization parameters of ^ ~ 10°-10^, and velocities 
from ^ 0-2000 km s~^. Our analysis of NGC 6860 indicates two components with ionized 
columns of ~ 10^°-10^^ cm~^, with ionization parameters of ^ ^45 and 180 ergs s~^, and 
outfiowing velocities between and 300kms~^. However, the constraints on the outfiow 
velocities are uncertain since measurements of individual lines uncover different central en- 
ergies depending on the model used. Thus, UV data will be particularly important to obtain 
for NGC 6860 to determine the kinematic properties of the outflowing material. 

We flnd complexity not only in the soft spectrum but also in the hard X-ray spectrum 
of NGC 6860. The hard spectra show a broad 11000 km s'^) FeKa line with EW ^ 100 
in both the XMM-Newton EPIC and Suzaku XIS observations. The FWHM of the i ron li ne 
is much higher than the measured FWHM for H/3 (5920±600kms~^: iBennert et al.l ( 120061 )). 



indicating that the Fe Ka line originates in a region interior to the broad line region and 
closer to the central black hole. The energy of this line is low (6.35 ± 0.02 keV), suggesting 
a red-shift in the 6.41 keV Fe I Ka line corresponding to f ~ 2800 km s~^. This energy is 
unusual, considering that typical measurements from Chandra high energy grating sp ectra 
have mean line-center energies of 6.399 ± 0.003 keV ( Yaqoob fc PadmanabhanI |2004| ). An 
additional possibly redshifted narrow Fe K emission line is signiflcant in the XMM-Newton 
spectra with an energy of 5.95lo;QgkeV (v ~ 0.08c) and EW = 17 eV. However, with the 
more physical diskline model, the 5.95 keV line is no longer signiflcant. 



The shape of the broad band continuum - measured from a combination of XMM- 
Newton, Suzaku, and Swift BAT observations - is well flt by an absorbed power law model. 
The column density of the neutral absorber is ~ 10^^ cm~^. The power law ind ex (P = 1.64) 



is lowe r than the typical value of 1.75 measured for the Swift BAT AGN sample (IWinter et al. 



2009aJ), but much higher than the value of 1.0 measured in the earlier short XMM-Newton 
observation of this source. The high energy cutoff for the power law is about 100 keV, which 
is consistent with a Comptonization spectrum. 
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Table 1. Features Detected in the Soft XMM-Newton EPIC spectra 



Component Energy (keV) 



a (keV) / (photons cm ^s ^) 



O VII edge 
O VIII edge 
Mg XI edge 
N VI hne 
N VII hne (?) 
Ne IX hne 
Si XIII hne 



0.719 

0.845: 
1.717 
0.405: 

n 4 p- 1+0.006 

Q qQ7+0.007 

-"-•o^^-O.Oll 



h0.004 
-0.004 
h0.004 
-0.003 
h0.039 
-0.023 
hO.006 
-0.405 



n 1 7C+0.013 

U-J-'O-o.ooy 



-0.005 
-0.019 



0.254: 

u.uz,/ _o.oo8 



0.01 
Q_01+0.176 

0.01 

0.01 



-4.47 



+0.37 
-0.57 



X 



;+3.84 



2.45_Lo:26 ^ 



10-^ 
10-^ 



-4.53 



+0.44 
-0.69 
_1 07+0.25 
-'-•'J' -0.23 



X 



10-^ 
X 10-5 



Listed are additional features that significantly (|Ax^| > 15) improved the fit to the EPIC 
pn and MOS spectra. These features were added in addition to a power law absorbed by 
the Galactic column in the 0.3-3 keV band. For the MOS spectra, we excluded data 
< 0.5 keV, due to inconsistencies with the pn spectra. Negative intensities indicate 
absorption lines, while positive intensities indicate emission lines. The recorded energies 

are in the rest frame of the AGN. 



Table 2. Lines Detected in the RGS spectra 



Energy (eV) /(xlO "^photons cm ^s ^) ID |Ax 

502i« i.nt'oil 

54511 -0.77l°f 



9171? -1-561°:^^ 



N VI 


17.7 


N VII 


48.7 


VII 


82.1 


Ne IX 


163.5 



Listed are the fitted energy, intensity, probable identification, and change in upon 
adding a Gaussian to model the indicated line at the redshift of the source {z — 0.014884). 
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Table 3. Spectral Fits in the 0.3-3 keV band with Warm Absorption Models 



Component 


Suzaku XIS 


XMM PN 


XMM RGS 


Nh"" (cold) 




6 41+0 04 

"■^-'--0.08 


q 7q+0-66 


N^^ (warm 1) 
log^^ (warm 1) 
vturb" (warm 1) 


' ■^'-'-3.41 
1 qr,+0.02 

i.yo_o.o9 

co+SO 


1 9 c:9+0.83 
-'-■^■"-'^-0.90 

2.04l°:°t 
2511 


14.1611:1 

9 A 9+0.05 
^■^^-0.13 
QQ+32 
O'J-26 


N//'' (warm 2) 
log^^ (warm 2) 
yturb" (warm 2) 


7.721^:^^ 

^•"'-'-0.07 




0.591°:!^ 
1.641°:^? 

113"^ 


kT (keV) 

T c 
-'feT 

r 

xVdof 


n in2+°-°22 

U.iUZ_Q 
•J-^' -0.14 

1-741°:°^ 

64.46l°:«t 
938.9/875 


n 1 oi +0.002 

U.iO±_o.o03 

27+°-°^ 
"-••^'-o.oi 

1 rr+0.01 
-■-■^^-0.01 

64.881?:^^ 
727.0/534 


0.87+ 

-■-•^"-0.39 

2.40l°:^« 

60.3211°^ 

2444.4/1951 



^Column density of the neutral (ztbabs) or warm absorber 
(warmabs) in units of 10^°cm~^. 

'^Additional components of the warm absorber model in- 
clude the ionization parameter (log^, where ^ has units of 
ergss"'^) and the turbulent velocity (in units of kms~^). 

'^Intensity of the model component in units of lO^"' pho- 
tons cm^^s^^ at IkeV for either the blackbody component 
(kT) or the power law (F) component. 

+The indicated parameter is an upper limit. 
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Table 4. Broad Band Spectral Fits in the 0.3-195 keV Range 



Component Fitted Value 



Nh'' (cold) 




kT (keV) 


U.UO/_oo09 


T b 


n 74+0-22 
•J- '^-0.17 


r 


-■-•"^-0.02 


Ecutoff (keV) 


103+3^? 


T b 


51.101°:^^ 


^FeKa (keV) 


a qq 1+0.030 


(JpeKa (keV) 


o.i2oi°:°i 


EW^eXa (eV) 


86.8lit^ 


T b 

i-FeKa 


n 00+0.04 

'-'■^^-0.04 


xVdof 


4254.3/3873 



Fo.3-iofcey (pn) 3.4 x 10"" ergssem-^ 
Fis-sofcev (PIN) 3.3 X 10-11 ergs s cm'^ 
Fi4-i95feey (BAT) 4.2 X 10" ergs s cm'^ 



'^Column density of the neutral absorber 
from the ztbabs model in units of 10^° cm"^. 

I' Intensity of the model component in units 
of 10"^ photons cm~2 at 1 keV for either 
the blackbody component (kT), the power 
law (r) component, or a Gaussian compo- 
nent. 
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Table 5. RGS Best-fit Absorption Values 



Parameter 


Fitted Value 


Nh^ (cold) 


A 24+0.90 


Nh^ (warm 1) 
log^^ (warm 1) 

Zabsorber (warm 1) 


2 26+0-24 
0:014269l°:S 


Nh'' (warm 2) 
logi^^ (warm 2) 

^absorber (warm 2) 


n fiq+O-17 

"-"•^^-0.12 

^-•"'-'-0.13 

n m 42^^1+0.000449 

U.Ui^ZOi_Q QQ0335 


xVdof 


2458.8/1952 



'^Column densities of the neutral 
and warm absorbers from the ztbabs 
or warmabs model are in units of 
10^0 cm-2. 

^The ionization parameter (log^, 
where ^ has units of crgss"^) is given 
from the warm absorption model. 
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Fig. 1. — The XMM- Newton pn background subtracted light curves are shown binned by 
1024 s (top) and 5760 s (bottom, normahzed by the average count rate in the specified energy 
band). Light curves are shown for three selected energy bands: 0.1-3keV (black), 3-7 keV 
(blue), and 7-10 keV (green). Clearly, the greatest amount of variability is seen in the soft 
band, which varies by about 2ctss~^. However, the normalized light curves (bottom) show 
that the relative amount of variability is the same in each of the bands. Further, each of the 
normalized light curves vary on the same timescale. We have excluded from the plot regions 
of high background fiares. 
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Fig. 2.— The Suzaku XIS (XIS0+XIS1+XIS3) normalized background subtracted light 
curves are shown in the soft band (0.1-3 keV, top), medium band (3-7 keV, middle), and 
hard band (7-12 keV, bottom). Each of the light curves are binned by 128s. As in the 
XMM- Newton observation, we find that variability in each of the bands is well-correlated. 
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Fig. 3. — The ratio of the data/model for the XMM- Newton EPIC spectra, fit by a simple 
absorbed power law from 3-10 keV, is shown for an absorbed power law model. The pn 
(black), MOSl (blue), and M0S2 (green) data are shown specifically around the Fe Ka 
emission line region. Clearly, the Fe Ka emission line is present. A narrow emission line at 
5.95 keV is also significant (Ax^ = —15). 
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Fig. 4. — Ratios of the data to model are plotted for spectral fits to the Suzaku observations 
in the 0.3-3 keV band. In the top plot, the ratio of the data/model for a simple absorbed 
power law model is shown. The fit is greatly improved with the addition of O VII and O VIII 
edges, as shown in the bottom plot. The Suzaku XIS0+XIS3 spectrum is shown in black 
and the XISl spectrum is shown in blue. Also present in the spectrum is an absorption edge 
at 1.91 keV, which is likely an instrumental feature associated with the Si K edge. 
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Fig. 5. — Ratios of the data to model are plotted for spectral fits to the XMM-Newton EPIC 
observations in the 0.3-3 keV band (rebinned by 50 for illustrative purposes). In the top plot, 
the ratio of data/model points is shown for a simple absorbed power law model. The fit is 
greatly improved with the addition of O VII and O VIII edges, in addition to additional lines 
described in Tabled! as shown in the bottom plot. 
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Fig. 6. — The ratio of the data/model for the power law + blackbody model, along with 
the O VII and O VIII edges, is plotted for the RGS spectra (re-binned to a signal-to-noise 
of 7 for illustrative purposes). The observed energy is plotted, while lines and features are 
labeled in the AGN rest-frame energy. The value where the data and model match (1.0) 
is plotted with a thick line. RGSl is shown in black and RGS2 is shown in blue. Regions 
where the RGSl or RGS2 have no data, are represented by a solid line. Emission/absorption 
features that were added to improve the spectrum (i.e. Ne IX near 0.9 keV, see Table [2]) 
are indicated. Additionally, there appears to be an Fe M-shell UTA feature in the spectra 
between 0.7 and 0.8 keV. Mg XI also appears to be present in absorption, though fitting a 
Gaussian at this energy does not significantly improve x^- 
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Fig. 7. — The ratio of the data/model for the two component warm absorber model described 
in § 14.2.41 are shown for the RGS spectra, re-binned to a signal-to-noise of 7. The observed 
energy is plotted, while lines and features are labeled in the AGN rest-frame energy. Overall, 
the two component warm absorber fit (along with a soft excess and power law continuum) is 
a good fit to the RGS spectra. Some of the worst fits to the data are at the lowest energies, 
particularly below ^ 0.7 keV. The region near O VII show that this region is the least well 
fit by the model. Similar results are seen in the XMM-Newton EPIC pn spectra. 
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Fig. 8. — Plotted is the unfolded spectrum for the broadband XMM- Newton pn, Suzaku XIS 
and PIN, and Swift BAT spectral fits. Our best fit model includes edges and Gaussian lines 
to fit the prominent warm absorption features (O VII, O VIII, N VII, and Ne IX), a black 
body to fit the soft excess, a neutral column of Nr ~ lO^^cm"^, a broad (a = 0.12 keV) Fe 
Ka line, and a cutoff power law. Best-fit parameters for this model are shown in Table HI 
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Fig. 9. — Plotted is the ratio of the data/model of the Fe Ka region using the Gaussian 
model (top panel) and a diskline model (bottom panel). Both the XMM-Newton pn (black) 
and Suzaku XIS0+XIS3 (red) spectra are shown. There is very little difference in the ratio 
plots from both models. Statistically, the reduced is slightly lower with the Gaussian 
model. 



